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Abstract: Scaling of mass ratios in intermediate volumes, obtained with im- 
proved SU(2) lattice actions is tested against analytic results for the Wilson 
and continuum action. A new improved action is introduced by adding a 
2x2 plaquette to the Symanzik action. Completing a square leads to a 
^ ! covariant propagator that simplifies perturbative calculations. Data is pre- 

sented on lattices of size 4 3 x 128, with lattice spacings of approximately 
0.02 and 0.12 fermi. For the latter case no further improvement as compared 
to the tree-level action was observed when including the Lepage-Mackenzie 
tadpole correction to the one-loop improved Liischer-Weisz Symanzik action. 
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O-t. 1 Introduction 

Improvement of lattice actions aims at doing Monte Carlo simulations on coarser lattices, 
such that with a modest number of lattice spacings the physical volume is sufficiently large. 
But perhaps more importantly it should make extrapolations to the continuum limit more 
reliable, as has been one of the main objectives in the non-perturbative determination of 
the running coupling constant [|TJ. Here we consider the Symanzik improvement scheme 0, 
which is designed to remove lattice artefacts by adding irrelevant operators to the lattice 
action, whose coefficients are tuned by requiring spectral quantities to be improved to the 
relevant order (on-shell improvement [Q]). For Symanzik improvement to work it seemed 
that unreasonably small values of the bare coupling constant were required. 

Mean field inspired Symanzik improvement || was introduced to beat the bad conver- 
gence of perturbative expansions in the bare coupling constant. In particular the Parisi 
mean field coupling || defined in terms of the plaquette expectation value is seen to im- 
prove considerably the approach to asymptotic scaling. Despite some attempts no good 
theoretical understanding for this is available. The tadpole prescription also includes cor- 
rections ]|] to the coefficients in the Symanzik improved action, which can be seen as a 
mean field renormalization of the link variables on the lattice. Only phenomenological ar- 
guments have been provided to support this. Standard tests in pure gauge theories involve 
restoration of rotational invariance H. More involved, but of direct physical relevance, are 
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the tests in charmonium spectroscopy @, used to extract a value of the strong coupling 
constant ||. 

We here stress the necessity of improving scaling, rather than asymptotic scaling, which 
in spectroscopy is less important since one has to set the scale by one of the masses or the 
string tension. In this sense our study is complementary to that of ref. 0. Although one is 
ultimately interested in the infinite volume limit, from the point of studying the approach 
to the continuum limit a finite volume provides a useful tool. If improvement fails there, 
it sheds doubt on results in large volumes (when successful, however, one does not imply 
the other). Perhaps a somewhat inappropriate comparison is that we consider our study 
as a well controlled laboratory experiment, where conditions are manipulated so as to rule 
out as much as possible external disturbances. 

The setup of this letter is to first introduce and motivate the new improved action. 
It simplifies certain perturbative calculations, and provides in part an analytic test of 
improvement in a small volume for which we present the effective potential in a constant 
abelian background field. Also the Lambda parameter of the new improved action is 
related to that of the Wilson action. We then present our Monte Carlo data at very small 
and intermediate volumes and end with conclusions. Details of the analytic study will be 
presented elsewhere (preliminary material and some further discussion can also be found 
in two communications to conferences |TTJ ) . 



2 Square Symanzik action 



There is a large redundancy in choosing an improved action, when parametrized in terms of 
Wilson loops. We shall use this to allow for a simplified "covariant" gauge choice, achieved 
by adding to the Liischer-Weisz (LW) Symanzik action a 2 x 2 plaquette. 
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The <> imply averaging over the two opposite directions for each of the links, called 



"clover" averaging in ref. [fLl}| . Numerical factors were chosen to agree with earlier conven- 
tions ||. Note that sometimes c<i and C3 are interchanged in the literature [ IT] , Here 
C4 is assigned to the 2x2 plaquette. 

The number of parameters required to improve the action to a certain order is simply 
determined from the number of gauge and hypercubic invariant operators that one can 
write down up to that order (read off from the dimension of the operator). For pure 
gauge theories there is only one operator of dimension zero and three of dimension two. 
One of these is redundant as it can be removed by a field redefinition, which can also be 
implemented at the level of the Wilson loop representation. It allows one to choose || 
c 3 = 0. 

As usual we relate lattice and continuum fields by U^{x) = Pexp(/ ° A^{x + sfi)ds). 
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This gives the following expansion for the lattice action [[El, [12| 
S({ci}) = - y(c + 8 Cl + 8c 2 + 16c 3 + 16c 4 ) E Tr 

X,fJ,,V 

+a 6 (| + c 3 ) E Tr(^F M (a;)^ A (x)) + a 6 | E Tr ((D^ A (x)) 2 ) 

x,fi,u,X x,/j,,u,\ 

+ ^(c + 20 Cl - 4c 2 + 4c 3 + 64c 4 ) E Tr (X^(:r)) 2 + 0(a 8 ) . (2) 

To fix the definition of the coupling constant one imposes (co + 8ci + 8C2 + 16c 3 + I6C4) = 1. 
Computing two particular spectral quantities as a function of these parameters allows one 
to determine these coefficients. At tree-level the conventional choice amounts to putting 
c = 5/3, C\ = —1/12 and c 2 = c 3 = c 4 = 0. The one-loop (0(g 2 )) correction to these 
coefficients was computed by Liischer and Weisz [[|. For c 4 7^ a similar calculation is 
in the process of being completed by one of us. At tree-level we have fixed c 4 by the 
following requirement. When expanding the action to quadratic order in the lattice field 
qn(x), defined by U^{x) = exp(q /M (x)), one finds 

$2 = E "I Tr [coiPtfM - duq^(x)) 2 + 2 Cl {(2 + d^d^x) - d v q^x))} 2 

X,fl,U 

+c 4 {(2 + d v ){2 + d^id^ix) - d u q^x))} 2 ] , (3) 
where <9 M is the lattice difference operator d fJj ip(x) = <p(x + fi) — <p(x). If we now choose 

c 4 = 2 2 c , z = d/c , (4) 
we can complete squares and obtain a simple gauge fixing function 

F gf {x) = v^E d l (l + *(2 + 5j)(2 + d,)) q,(x). (5) 

It is for this reason we propose to call the new improved action the square Symanzik action. 
At tree-level one finds 

c = 16/9, c a = -1/9, c 2 = 0, c 3 = 0, c 4 = 1/144. (6) 

An amusing, and potentially useful feature is that the relation c 4 Co — c 2 is not affected by 
tadpole corrections, where one replaces U^{x) by U^(x)/u , with u the fourth root of the 
average value of the plaquette. 



u 



^Re^Tr (Q). (7) 

Here k is the number of colors. For values of uq 7^ 1 one easily finds z = — 1/(16wq) and 
c = l/(l + 4z) 2 . 

In the covariant gauge the propagators for the ghost and vector fields are simply given 

by 

Ghost : P(k) -- 



v/co Ea (4 sin 2 (A; A /2) + Az sin 2 k } 
Vector : P^vik) = 



P(k)8^ 



co~(l + 4zcos 2 (^/2)) • 
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It illustrates an important feature of improved actions with more than nearest neighbor 
couplings in the time direction: unphysical poles appear at the scale of the cutoff, 1/a. For 
low energy physics they are harmless in the same way Pauli-Villars regulator fields are 
harmless at energies below the scale of the cutoff. However, on the lattice these spurious 
poles are more cumbersome to handle as they do not simply appear in loops (i.e. vertices 
do not preserve something like ghost number). Nevertheless, there is a way of separating 
off their contributions [ 10] . Each of the propagators can be factorized in the sum of two or 



three standard (single pole) lattice propagators, P s = 1/(4 sin (§/co) +WsW), 

P(k) = Z(k)(P-(k)-P + (k)), (9) 

p, v {k) = s, v (z;(k)p^k)-z+(k)p + (k) + z^k)p (k)), 

It is straightforward to derive the explicit expressions for the Z factors and energies u from 
eq. (8). Note that Z° = and that (for u = 1) cj 2 (0) = uo 2 _ (0) = 12. The spurious poles 
in this case occur at an energy 2asinh(v / 3)/a. 

One particularly simple test of improvement is achieved by computing for SU(2) the ef- 
fective potential for a static abelian zero-momentum background field, Uj = exp^iCja^/N) 
and Uq = 1, that is a solution of the (lattice) equations of motion. Introducing the quan- 
tum fluctuations through U^(x) = e^w^ one easily diagonalizes the quadratic fluctuation 
operator in the covariant gauge 

? B f = v^E t>l + *(2 + Dl)(2 + D,)) (10) 

where D^(x) = U^ip{x + jl)W — (p(x). Due to the background field, momenta will be 
shifted to k s = (2irn + sC)/N, where s = for the isospin neutral and s = ±1 for the 
isospin charged components of quantum fields. The eigenvalues can be directly read off 
from eq. (8) and one finds 



V f°{C) = N E log + 4asinh 2 «c 



u 2 




l+4z + — + u\ 1 + — , 11 



with Xj = l+4zcos 2 (iA;+) and u 2 = J2j 4A^ sin 2 (ifc+). The result, normalized to Vf 3 {0) = 0, 
is plotted in figure 1 for u = 1 (together with the effective potential for the Wilson action, 
obtained by taking z = 0). At iV = 6 we can not distinguish the result from the continuum 
at the scale of this figure. 

We can use the abelian background field also to compute the one-loop correction to 
the tree-level kinetic term ±gQ 2 (dCi(t)/dt) 2 , which yields \{g~ 2 + ai)(dCi(t)/dt) 2 . In the 
continuum limit g~ 2 = g^ 2 — 11 log(iV)/127r 2 is kept fixed while sending the number of 
lattice spacings, N, to infinity. An analytic expression for a\(N) was found in terms of 
a sum over spatial momenta, which reduces to the result for the Wilson action at 
2 = 0. Computing this sum for one hundred lattices and fitting the result to a polynomial 
in l/N we find a x = -0.0340012235(1) for z = -1/16 and a x = -0.1648688946(1) for 
z = 0. From the difference one determines the ratios of the Lambda parameters between 
the square Symanzik (A52) and the Wilson (A w ) actions. One can also compute the 



4 



one-loop correction for Jg^Tri^-, which gives an identical result for the Lambda ratios. 
Alternatively, we used the heavy quark potential method [1^], which also allowed us to 
extract the Lambda ratios for SU(3). We quote the following result 



A S 2/A W [SU(2)] = 4.0919901(1), A 52 /A H/ [SU(3)] = 5.2089503(1). (12) 

In addition, the one-loop perturbative expansions of the SU(fc) expectation values for an 
a x b plaquette, U(P(a, b)), are given by 

(Re k~ l Tr U(P(a, &))) = 1 - \gl{k - *T b )- ( 13 ) 

For the square Symanzik action (eq. (6)) we find 

1) = 0.3587838551(1), w(l, 2) = 0.6542934512(1), w(2, 2) = 1.0887235337(1). 

(14) 



3 Monte Carlo data 

We wish to determine in small volumes the mass for the scalar (Af) and tensor glueballs, 
the latter split due to the breaking of rotational invariance in the doublet E + and the 
triplet T 2 + . Also the energies of the electric flux ("torelon") states with one, two and 
three units of electric flux (e^, i = 1,2,3) will be measured. In addition we consider the 
states with two (T-q or 5(110)) and three (T 2 + (lll)) units of electric flux that have T 2 + 
quantum numbers (negative parity in two directions of electric flux, symmetrized in those 



two directions). See ref. [0| for details and further references. 

The size of the lattice used is 4 3 x 128 and masses m are converted to dimensionless 
parameters into z = mL; in lattice units we hence multiply the mass with the number of 
lattice sites in the spatial directions. In large volumes one should have z ek = aL 2 \/k } where 
a is the infinite volume string tension. This is why we will consider the rations Jz ek / z A + . 



These and other mass ratios will be plotted as a function of z A +. The analytic result fTJ 



derived by diagonalizing an effective Hamiltonian to describe low-lying states is valid up 
to z A + ~ 5, after which degrees of freedom that were integrated out perturbatively will 



receive non-perturbative contributions The breakdown will occur at smaller volumes 
for higher excited states. 

For the Wilson action we have chosen (3 = 3.0 and (3 = 2.4; for the improved actions 
(3 was tuned to yield results in roughly the same physical volume. These parameters 
correspond to lattice spacings of approximately a = 0.018 and a = 0.12 fermi. For the 
smallest of these two, one expects tree-level improvement to be effective and we have 
therefore not tadpole corrected the actions in this case. Note that for these small volumes 
one finds from the analytic results that the lattice artefacts in the mass ratios are quite 



much bigger [TJj than in larger volumes. Data was taken for both the LW and square 



Symanzik actions, and as a test on our programs also for the Wilson action for which we 
can compare with available high precision data ||15|| . 

At the larger volume we concentrated our attention to the LW Symanzik action with 
tree-level and tadpole corrected one-loop values of the coefficients. We verified that there 
is no observable volume dependence of uq by comparing its value with the one on an 8 3 x 64 
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lattice (the difference was less than 0.3%, consistent with zero within statistical errors). 
Following the prescription of refs. ||, [?|, Q we took for SU(2) 

c = 5/3, ci = -(1 + 0.2227a s {u ))/{12u 2 ), c 2 = -0.02224 x 5a a /(3ug), c 4 = 0, 

a s (u ) = -(41ogu )/1.725969, (15) 

obtained from the one-loop coefficients determined by Liischer, Weisz and Wohlert |3|, 
12| . Substituting these coefficients in eq. (1) and multiplying by /3/4 = 1/g 2 . gives the 
action we used for our simulations. We do not absorb the tree-level value of Co in the 
definition of the coupling constant, as was done in ref. |7], [|. When using the convention 
of eq. (15) the standard two-loop relation between (3 and aA needs no modification. But 
the Lambda parameter has to be corrected for the fact that the Liischer- Weisz choice of 
coupling amounts to multiplying eq. (1) by (g^ 2 + 0.08112), so as to compensate for the 
one-loop correction to Co. 

In the intermediate volumes we have used for our simulations, masses remain small com- 



pared to the spurious unphysical poles. This allows us to use the variational approach [16 
to increase the overlap of the the states to be measured. We have been able to extract 
clean signals. On very coarse lattices where masses would no longer be small in lattice 
units, one looses the signal in the noise too early to extract it reliably, whereas also the 
variational method is no longer well founded. Recently these problems were tackled by 



using anisotropic lattices ||, well known from finite temperature studies \T%. Only imple 



menting improvement for the spatial directions will in addition remove the problems with 
a non-hermitian transfer matrix. 

The raw data are listed in tables 1 and 2, based on performing the variational analysis 
on the second time slice. We have verified that the result is stable against performing the 



variational analysis on the first time slice. We used 3 to 8 operators, as defined in ref. [15] 



for the variational analysis. They were computed in terms of Teper-fuzzed links QTHJ . Only 
for the determination of the scalar glueball mass at a ~ 0.02 fermi the variational analysis 
was important, in most other cases a single but Teper-fuzzed operator was sufficient to 
obtain accurate results. 

Another issue is that for a ~ 0.02 fermi the small value of the coupling gives rise to 
large autocorrelations that can affect the energies of electric flux. In most cases we found it 
useful to correct for this by eliminating data for which the average of the spatial Polyakov 
loops over the 128 time slices (and a few heat bath updates) was bigger in absolute value 
than one half. Our results for the Wilson case at /3 — 3 agree to high accuracy with those 
reported by Michael |T5| . 

Because of the availability of analytic results, it is not necessary to exactly tune the 
different actions to the same physical volume. Nevertheless in particular for the data at 
a = 0.12 fermi we made an effort to tune parameters appropriately, as we can make a 
stronger point when directly comparing lattice data at the same physical volume. The 
value of Mq is determined self consistently ]7|, [5], adjusting with the help of the Ferrenberg- 



Swendsen trick |19| the input value of u to agree with its measured value. This only 
requires little Monte Carlo time. The results of tables 1 and 2 are presented in figure 2 
to compare with the analytic results for the continuum (solid curves) and for the Wilson 
action on a lattice of size 4 3 x oo. We have used approximately 160 hours of CPU time 
on a Cray C98 to generate and analyze the data presented in this paper. Computational 
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overheads for improved actions amount to a factor 3 for the LW and 4 for the square 
Symanzik action over the standard Wilson action. 

4 Discussion 

As was to be expected, at lattice spacings of approximately 0.02 fermi (z A + ~ 2), tree- level 
improvement is seen to bring the lattice results quite close to those of the continuum, both 
for the LW and square Symanzik actions. In both cases the improvement is considerable. 

Also at lattice spacings around 0.12 fermi and volumes of approximately 0.48 fermi 
(z A + ~ 4), the agreement of the Wilson action lattice data with the corresponding analytic 
results is in general very good for the lowest lying states. The difference in the analytic 
result between the continuum and Wilson lattice action gives an indication how far the 
improved data is removed from the continuum result. Significant improvement is observed 
in some of the cases, in particular for z T + /z A +, approaching the continuum analytic result. 

The most salient feature of our data is that tadpole correction has no significant effect 
on the tree-level improved data for the ratios. Perhaps for the cases where tree-level 
improvement is already significant this is what one would want, but our results show some 
instances where tree-level improvement has no effect and the tadpole correction is of no 
help either. 

In particular we note that the ratio yfz^j z A + , measured to an accuracy of better than 
1.5%, deviates from its continuum value by 5-6%. For this quantity tree-level improvement 
as well as tadpole corrected one-loop improvement are unable to show deviations from the 
Wilson result. This result puts some doubt on the usefulness of the tadpole correction for 
careful extrapolations of mass ratios to the continuum limit. 

One might object that the lattice spacing we have used to implement the tadpole 
correction, a = 0.12 fermi, is not really large enough. We have certainly not probed lattice 
spacings as large as a = 0.4 fermi, that have been advertised 0. Nevertheless for a = 0.12 
fermi, Uq = 0.6819(1) and significantly deviates from 1. The correction to C\ at these 
parameters is 27% with respect to its tree-level value (without tadpole correction it would 
have been 17%). 
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Wilson action, (3 = 3.0 



Rep. 


#op- 


1/0 


2/1 


3/2 


4/3 


At 


5 


2.105(9) 


2.03(1) 


2.02(2) 


2.08(3) 


E + 


3 


1.743(5) 


1.703(9) 


1.71(1) 


1.70(2) 




3 


3.315(9) 


3.25(2) 


3.21(4) 


3.2(1) 


ei 


3 


0.277(3) 


0.269(5) 


0.269(6) 


0.270(7) 


e 2 


3 


0.588(5) 


0.575(7) 


0.576(9) 


0.58(1) 


e 3 


3 


0.978(5) 


0.962(8) 


0.97(1) 


0.97(2) 


LW Symanzik action, (3 = 2.374 


Rep. 


#op. 


1/0 


2/1 


3/2 


4/3 




6 


1.89(1) 


2.01(2) 


2.03(3) 


2.09(5) 


E + 


3 


1.626(6) 


1.77(1) 


1.78(2) 


1.80(3) 


J 2 


3 


3.134(8) 


3.43(2) 


3.47(5) 


3.4(1) 


ei 


3 


0.307(3) 


0.334(5) 


0.339(6) 


0.342(8) 


e 2 


3 


0.656(5) 


0.718(8) 


0.73(1) 


0.74(2) 


e 3 


3 


1.104(6) 


1.21(1) 


1.24(1) 


1.26(2) 


Square Symanzik action, /3 = 2.2013 


Rep. 


#op. 


1/0 


2/1 


3/2 


4/3 


At 


6 


2.15(1) 


2.31(2) 


2.30(3) 


2.29(5) 


E+ 


3 


1.844(6) 


2.00(1) 


2.02(2) 


2.04(3) 




3 


3.56(1) 


3.87(3) 


3.88(7) 


3.8(3) 


ei 


3 


0.375(4) 


0.405(6) 


0.408(7) 


0.412(9) 


e 2 


3 


0.808(6) 


0.877(9) 


0.89(1) 


0.90(2) 


e 3 


3 


1.368(8) 


1.49(1) 


1.52(2) 


1.54(9) 



Table 1: Values of z — mL at a lattice spacing of approximately 0.02 fermi for SU(2) on 
a 4 3 x 128 lattice. We have performed 16000 measurements (25 heat-bath sweeps apart) 
for Wilson and 20000 for both LW and square Symanzik actions (10 sweeps apart). The 
entries in the table correspond to the representations of the cubic group, the number of 
operators used in the variational analysis and the effective masses extracted from n/l ratios 
of correlation functions, i.e. - log (£t < 0(t + n)0{t) > / £t < 0(t + £)G(t) >). Entries 
in boldface are taken as final estimates for figure 2. Errors have been analyzed using the 
jackknife method. 
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LW Symanzik action, (3 = 1.83 



Rep. 


#op. 


1/0 


2/1 


3/2 


4/3 


At 


7 


3.71(1 ) 


3.74(2) 


3.78(5) 


3.9(2) 


E+ 


7 


3.212(9) 


3.29(2) 


3.30(4) 


3.3(1) 


T+ 


3 


6.13(1) 


6.31(6) 


6.3(3) 


6.2(1.0) 


ei 


7 


0.813(6) 


0.84(1) 


0.84(1) 


0.84(2) 


e 2 


7 


1.75(1) 


1.80(2) 


1.80(3) 


1.80(3) 


e 3 


8 


2.89(2) 


3.05(3) 


3.09(5) 


3.2(1) 




7 


1.857(6) 


1.92(1) 


1.92(2) 


1.91(3) 


T+(lll) 


7 


2.67(1) 


2.67(2) 


2.66(3) 


2.67(7) 


Tadpole corrected LW Symanzik action, j3 = 


2.04 


Rep. 


#op. 


1/0 


2/1 


3/2 


4/3 


At 


7 


4.07(1) 


4.06(3) 


4.05(8) 


3.8(2) 


E+ 


7 


3.366(6) 


3.57(2) 


3.57(5) 


3.6(1) 




3 


6.28(2) 


6.76(7) 


6.6(4) 


6.5(1.0) 


ei 


7 


0.889(5) 


0.94(1) 


0.94(1) 


0.94(2) 


e 2 


7 


1.920(1) 


2.02(2) 


2.01(3) 


1.98(4) 


e 3 


7 


3.402(14) 


3.41(3) 


3.38(4) 


3.21(7) 


Tit 


7 


1.893(6) 


2.06(1) 


2.07(3) 


2.07(4) 


T+(lll) 


7 


2.79(1) 


2.87(2) 


2.86(4) 


2.83(8) 



Table 2: The same as in table 1 but for a lattice spacing of approximately 0.12 fermi. We 
have performed 40000 and 48000 measurements respectively for tree-level LW and tadpole 
corrected one-loop improved LW Symanzik actions, in both cases separated 2 heat-bath 
sweeps apart. 
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V 




Figure 1: The SU(2) effective potential for a constant abelian background field C = 
(C, 0, 0). The full line represents the continuum result (obtained by taking the number of 
lattice spacings N — > oo). The lower two dashed curves are for the square Symanzik action 
(uo = 1) with N = 3 and 4. The upper three dotted curves are for the Wilson action with 
N = 3,4 and 6. 



11 



0.6 




0.3 



0.2 - 




1.8 



1.7 



1.6 



1.5 



*2 H l 



J I I I I I I I I I I I I L_ 





0.6 



0.5 



\ z tj(iii)/ z a; 




2t+ / z. t 
1 it rt i 



j i i i i i i i i i i i i i_ 



Figure 2: SU(2) Monte Carlo data (see tables) on a 4 3 x 128 lattice for the Wilson action 
(circles for our data and crosses for data by Michael [15], with tilted error bars when 
data overlap), the LW Symanzik improved action (triangles), the square Symanzik action 
(squares) and the tadpole corrected one-loop LW Symanzik action (pentagons) at lattice 
spacings of approximately 0.02 and 0.12 fermi. A comparison is made with analytic results 
for the continuum (solid lines) and Wilson action on a lattice of size 4 3 x oo (dashed lines). 
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